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REALITY PROBLEMS IN THE SOLITON THEORY 



P.G.GRINEVICH AND S.P.NOVIKOV 

Abstract. This is a survey article dedicated mostly to the the- 
ory of real regular "finite-gap" ( algebro-geometrical) periodic and 
quasiperiodic Sine-Gordon solutions. Long period this theory re- 
mained unfinished and ineffective, and by that reason practically 
had no applications. Even for such simple physical quantity as 
"Topological Charge" no formulas existed expressing it through 
the "Inverse Spectral Data". Few years ago the present authors 
solved this problem and made this theory effective. This article 
contains description of the history and recent achievements. It 
describes also the reality problems for several other fundamental 
soliton systems. 



1. Introduction 

The most powerful method for constructing exphcit periodic and 
quasiperiodic solutions of soliton equations is based on the Finite- 
gap or Algebro-geometrical approach, developed by Novikov (1974), 
Dubrovin, Matveev, Its, Lax, McKean, Van-Moerbeke (1975) for 1 -|- 1 
systems, extended by Krichever (1976) for 2 -|- 1 systems like KP (see 
EH, Eni, 0, El, Ei for extra information). Already in 1976 new 
ideas were formulated how to extend this approach to the 2-1-1 sys- 
tems associated with spectral theory of the 2D Schrodinger operator 
restricted to one energy level (see EHl )• These ideas were devel- 
oped in 1980s by several people in the Moscow Novikov's Seminar (see 
below). The "Spectral Data" characterizing the associated Lax- type 
operators consist of Riemann Surface ( "Spectral Curve" ) equipped by 
the selected set of points ("Divisor of Poles", "Infinities"). In the fi- 
nite gap case this Riemann surface has finite genus, and the number 
of selected point is also finite. The algebro-geometric approach in par- 
ticular allows to write down explicit solutions in terms of the Riemann 
6'-functions. 

In modern literature very often the problem is assumed "more or 
less" completely solved if such formulas are derived. However, in some 
cases this belief is too naive and does not correspond to the needs of real 
life. For example, it is necessary to select physically or geometrically 
relevant classes of solutions corresponding to the source problem (i.e. 
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solutions, satisfying some reality conditions, regular solutions, bounded 
solutions and so on). Is it easy or not? 

To reach this goal, following problems should be solved. 

• Problem 1. How to select solutions real for real {x,t)7 

• Problem 2. How to select real nonsingular solutions? 

• Problem 3. How to select periodic solutions with given period 
(or quasiperiodic solutions with given group of quasiperiods)? 

Remark. We call solution non-singular, if it is non-singular on the 
whole real Abel torus. It should remain non-singular under action of 

all (real) higher flows from the corresponding integrable hierarchy. 

The generic x-direction is normally ergodic in the Abel torus, so this 
definition is equivalent to the standard one. However, for some spe- 
cific values of constants of motion theoretically we may have solutions, 
which are regular in the standard sense, but blow-up under the action 
of the higher symmetries. 

For some models like Korteveg-de Vries equation (KdV), defocus- 
ing Noninear Schrodinger Equation (NLS), Kadomtsev-Petviashvity 2 
(KP2) equation selection of real and nonsingular sohition is rather 
straightforward. But for many other models like focusing NLS, Sine- 
Gordon equation (SG), KPl, inverse scattering transformation for 2-D 
Schrodinger operator based at one energy the problem of selecting real 
solutions is rather difficult. 

The theory of ^-functions is very complicated and ineffective. The 
complexity is hidden behind the simple notations in these formulas. 

Our goal is to discuss in more details the Sine- Gordon equation: 

(1) utt - Uxx + smu{x,t) = 0. 
In the light-cone coordinates 

(2) x = 2(^ + 77), t = 2(e-77), 
it has the form. 

(3) u^jj — Asinu, u — u{^,r]). 

According to our definition, the solution u{x,t) is x-pcriodic with pe- 
riod T if exp{iu} is x-periodic with that period. For the function u we 
have 

u{x + T,t) = u{x, t) + 27rn, n G Z. 

We call the quantity n a Topological Charge corresponding the the 
period T. 

We call the ratio n/T a, Density of Topological Charge. 
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The Density of Topological Charge can be naturally extended to all 
real generic regular finite-gap (quasiper iodic) solutions. It is the most 
basic conservation law. 

Problem: How to calculate topological charge of the real finite-gap 
solutions in terms of the spectral data. 

Let us remind, that the Inverse Scattering (Spectral) Data for KdV 
and Sine- Gordon systems consist of a Riemann surface (Spectral Curve) 
r with finite genus equals to g and a collection of points ("divisor") 
D = ji + . . . + jg. (For NLS and some other systems number of poles 
may be different from genus). 

In the case of KdV (or finite-gap periodic Schrodinger operator L = 

—d1 + u{x)) this surface V is hyperelhptic. In the case of the Sine- 

29 

Gordon equation the surface is also hyperelhptic, /i^ = A n(A — Aj) 

1=1 

with branching points (0, Ai, . . . A23, 00). However the classes of ad- 
missible Riemann surfaces and divisors for KdV and Sine- Gordon are 
dramatically different (see below). 

The 6'-functional formulas for Sine-Gordon were obtained by Kozel, 
Kotlyarov and Its [2215 ^H] in 1976. The reality problem remained 
unsolved. Indeed, the class of admissible Riemann surfaces was found in 
these works (see [I^). The nonzero finite branching points (Ai, . . . \2g) 
can be either real negative (Ai, . . . A2fc) G M or complex conjugate with 
nonzero imaginary part A2fc+i = A2fc+2, • • • , A23-1 = \2g- However, 
no ideas were proposed where the poles are located on the Riemann 
Surface. 

As it was realized in early 1980s by McKean Dubrovin, Novikov 
and Natanzon [H], [7j, Ercolani, Forest and McLaughlin PHIjCII! this 
problem is nontrivial. By that reason periodic finite-gap Sine-Gordon 
theory long period had no applications. 

An important idea how to describe position of poles for the real non- 
singular solutions was in fact suggested by Cherednik in 1980 j^. He 
was the first author who discovered (ineffectively) that for the given 
admissible real Riemann Surface there can be many different real Abel 
tori generating real nonsingular quasiperiodic solutions. Their number 
is equal to 2^ where 2k is the number of negative real branching points. 
All real finite-gap solutions are nonsingular for Sine-Gordon for the 
generic Riemann surface. His work was written in the abstract algebro- 
geometric form, and he never developed these ideas later. Extending 
this approach on the basis of "Algebro- Topological" ideas, Dubrovin 
and Novikov P| presented an interesting idea how to calculate topo- 
logical charge in terms of the "inverse spectral data". However, as it 
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was pointed out by Novikov in 1984 [S2], there was mistake in their 
arguments: the formula proposed in jU] was proved only for the small 
neighborhood of some very special solutions. The problem remained 
open till 2001. The complete solution (confirming Dubrovin- Novikov 
formula) was obtained by the authors in [13j as a development of the 
"Algebro-topological approach", suggested in P, see also [131, [ISl- 
is interesting that in the works j7j and later TD^, these components 
were described as the real subtori in the Jacobian variety J{T). How- 
ever this "^-functional description" did not led yet to any formula for 
the topological charge. It does not require any specific basis of cycles. 
As we know now, good formula for the topological charge can be writ- 
ten in very specific basis only. We believe that using this basis of cycles 
one can deduce our formula from the ^-functional expression. It would 
be good to do that. 

2. Physically relevant classes of solutions for the 
different soliton systems 

The Korteveg-de Vries equation (KdV) 

(4) ut + u^xx - Quux = 0, u = u{x, t), 

was originally derived in the water waves theory. As it was discovered 
in early 1960s (see introduction to the book [S3]), it naturally appears 
as a first non- vanishing correction for the dispersive nonlinear systems 
if dissipation can be neglected. In these models only real non-singular 
solution are physically relevant. 

Integration of KdV equation is based on the "Inverse Scattering 
Transform" for the 1-dimensional Schrodinger operator 

(5) L = ~dl + u{x,t). 

Selection of real KdV solutions is rather straightforward. 

(1) The spectral curve F /x^ = i?2g+i(A) should be real. It means, 

that -R2g+i(A) = A^^+^ + J2Pi^^ real coefficients Pi G M, or 

i=0 

equivalently, all roots are either real or form complex conjugate 
pairs. Therefore we have a holomorphic involution r : (A, n) 
(A, — /2) on r. 

(2) The divisor D should be real with respect to r: tD = D, or 
equivalently, the unordered set of points 7i,. . . , jg is invariant 
with respect to r. Of course, r may interchange some of them. 

Real nonsingular KdV solutions correspond to the following special 
spectral data: 
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(1) All branching points of Xk of F are real and distinct. Assume, 
that Ai < A2 < . . . < X2g+i- Then r has exactly g + I real 
ovals over the intervals Oq = {—00, Ai], Oi = [A2, A3], ... ,ag = 

(2) Each finite oval a^, 1 < k < g contains exactly one divisor point 

Ik e Ofc. 

Remark. A real curve of genus g may have at most g + 1 real oval. 
Curves with g + 1 real ovals (maximal possible number) are called M- 
curves. 

Generic finite-gap solutions are quasiperiodic with g incommensu- 
rable periods. How to select a;-periodic solutions with prescribed pe- 
riod T? Avoiding any use of Algebraic Geometry and Riemann Sur- 
faces, nice approach to the characterization of the strictly x-periodic 
solution in terms of the so-called "quasimomentum map" was devel- 
oped by Marchenko and Ostrovskii in 1975 |27j. This map was studied 
in the Quantum Solid State Physics Literature in 1959 (see It 
is well-defined in the upper half-plane outside of some vertical edges. 
Its analytical properties were effectively used in [77j. For example the 
approximation of x-periodic solution (potential) by the finite-gap ones 
periodic with the same period, was proved. Another approach based 
on the so-called isoperiodic defomations of finite-gap potentials was de- 
veloped by Grinevich and Schmidt in 1995 ^Tj. In the KdV case the 
isoperiodic deformations can be interpreted as the so-called "Loewner 
equations" for the corresponding conformal map. Let us point out, 
that there exists a big literature, dedicated to the KdV solutions with 
real poles (rational solutions, singular trigonometric and elliptic solu- 
tions) - see X where these ideas were started. These solutions are 
very important from the mathematical point of view: for example, the 
dynamics of poles satisfies to the equations of the rational and elliptic 
Moser-Calogero models respectively. However, they are related neither 
to nonlinear wave problems nor to the spectral theory of the corre- 
sponding Schrodinger operators. So we do not discuss this literature 
in the present survey article. 

The modified Korteveg-de Vries equation has the form: 

(6) vt + Vxxx - Qv'^Vr, = 0, V = v{x,t). 

It is connected with KdV by the Miura transformation: 

(7) u{x,t) = Vx{x,t) + v'^{x,t). 

The real non-singular solutions are physically relevant. 

The "complex" Non-linear Schrodinger equation (NLS) is a system 
of equations for the pair of independent complex functions q = q{x,t), 
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r{x, t): 

iqt + Qxx + 2gV = 



-irt + r^x + 2r^g = 

This system has 2 natural real reductions: defocusing NLS: r{x,y) 
-q{x,y) 

(9) iqt + Qxx -'^\q\'^q = 0, 



and self-focusing NLS r(a;, = y) 

(10) igt + g,, + 2|g|2g = 0. 

These equation describes nonlinear media with dispersion relations de- 
pending on the square of the wave amplitude (see [SS])- Among the 
todays applications of NLS is the theory of light propagation in the 
fiber optics. The sign + or — is determined by the dispersion relation, 
and the qualitative behaviour critically depends on it. From the math- 
ematical point of view, the defocusing NLS system is much simpler 
because the linear Lax operator is self-adjoint. The focusing NLS is 
much more complicated. In both cases physical applications requires 
regular solutions. 

The complex NLS spectral data are following: A hyperelliptic Rie- 
mann surface F with 2g + 2 finite branching points Ai,. . . ,X2g+2 and 
g + I divisor points D = 71 + . . . + 7g+i. In contrast with the KdV 
case, there is no branching at 00. 

Solutions of the defocusing NLS correspond to the following spectral 
data: 

23+2 

(1) F is real, i.e. the polynomial -R2C/+2 = 11 (-^ ^ ^k) has real 

k=l 

coefficients. F is defined by fi^ = R{X). The anti-holomorphic 
involution on F is defined by the map r : — where 
(A,/i) (A, -/i). 

(2) The divisor D is real with respect to r: tD = D. 

Selection of regular solutions is also very similar to the KdV case 

(1) All branching points of F are real. Therefore F has g + 1 real 
ovals over the intervals [A2fc-i, A2fc], k = 1, . . . , g + 1, i.e. F is an 
M-curve. 

(2) There is exactly one divisor point at each real oval. 

Selection of x-periodic solutions is completely analogous to the KdV 
case. 

Let us describe the data generating real solutions of the self-focusing 
NLS equations for regular spectral curves. By the Cherednik theorem 
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these solutions are automatically non-singular. Solutions, corre- 
sponding to singular spectral curves can be obtained as proper degen- 
erations. In contrast with the defocusing case, singular curves may 
generate regular x-quasiperiodic solutions. 

(1) r is a real hyperelliptic surface of genus g with 2(yf + 2 Ai < A2 < 
. . . < )^2g+2 finite branching points. There are no branching 
points on the real line, so they form complex conjugate pairs. 
The antiholomorphic involution r acts on the A-plane as rA = A. 
The points in F lying over the real line are invariant with respect 
to r. Equivalently, r : (A,/i) (A,/i). 

(2) There exists a meromorphic differential Vt such, that: 

• = (1 + o{l))d\ at the infinite points of F. 

• f2 is regular outside infinity. Therefore it has exactly 2g + 2 
zeroes. 

• Let D = 7i + . . . + 7g+i. Then the divisor of zeroes of is 
D + tD. Therefore, D + tD = 2ooi + 2oo2 - K. 

The Sine-Gordon equation in the light-cone variables was derived in 
the end of the 19-th century. It describes immersions of the negative 
curvature surfaces into M^. Assume, that an asymptotic coordinate 
system is chosen (a coordinate system such, that coordinate lines have 
zero normal curvature). The angle between the coordinate lines satisfy 
(jSI). It means, that only real regular solutions such that u{x,t) 7^ 0( 
mod tt), are relevant. 

The Sine-Gordon equation describes also dynamics of the Josephson 
junctions. In this model u{x,t) is the phase difference between the 
contacts, therefore the real non-singular solutions are relevant. How- 
ever, according to the leading experts in the Superconductivity Theory, 
the problem always requires boundary problem, so we have to consider 
either the finite interval or the half-line. 

The elliptic Sinh-Gordon equation 

(11) Uxx + Uyy + 4:H sinhu = 0. 

describes the constant mean curvature surfaces with genus equal to 
one, outside umbilic points (see review PJ). The constant mean cur- 
vature tori has no umbilic points, therefore real nonsigular solutions 
should be selected. In contrast with soliton equations, all real smooth 
double-periodic solutions are automatically finite-gap here: jUI -P^- 
It follows from the following observations by Hitchin JTH^: all isospec- 
tral fiows from the corresponding hierarchy are zero eigenf unctions of 
the linearized problem. But the linearized system is the 2-dimensional 
(elliptic) Schrodinger operator, and it may have only finite-dimensional 
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space of double-periodic zero eigenfunctions. It means, that the hier- 
archy contains only finite number of linearly independent flows at this 
point. As a corollary the spectral curve has finite genus. A further de- 
velopment of this idea was used by Novikov and Veselov in the paper 
|34j . It was shown, that all periodic chains of Laplace transformations 
consisting of the 2-dimensional double-periodic Schrodinger operators 
with regular coefficients are algebro-geometric (2D analogs of "finite- 
gap" operators). 

The Boussinesq equation 



is used for describing the water waves. For physical applications it 
is necessary to select real non-singular solutions. We would like to 
point out, that the problem of selecting such solutions in terms of the 
finite-gap data remains open. 

The Kadomtsev-Petviashvili (KP) equation 

(13) {ut + u^^^-6uU:„)x + 3a^Uyy = 0, u = u{x,y,t), G R. 

The auxiliary linear operator for KP has the form 



If a is imaginary, we have the so-called KPl equation, and L is the 
one-dimensional non-stationary Schrodinger operator. If a is real, L is 
the parabolic operator. In both cases the real non-singular solutions 
are physically relevant only. The necessary and sufficient conditions 
for the finite-gap spectral data selecting the real non-singular solutions 
were found by Dubrovin and Natanzon jH]. 

Real non-singular solutions of the KP-2 equation correspond to the 
following geometry: 

(1) r is a algebraic surface of genus g with a marked point and 
an antiholomorphic involution r such, that the marked point 
is invariant under the action of r. The marked point is the 
essential singularity of the wave function. 

(2) T has exactly g + 1 fixed oval, i.e. F in an M-curve with respect 
to r. Denote the oval containing the essential singularity by qq 
and the other ovals by a„, n = 1, . . . , g. 

(3) Each oval a„, n 7^ contains exactly one divisor point. 

In the case of the Kadomtsev-Petvishvili 1 equation the reality con- 
straints on the spectral curve are exactly the same as in the KP 2 case, 
but the divisor D has a completely different description: There exists a 
meromorphic differential Q with exactly one second-order pole located 





(14) 
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at the marked point such, that the divisor of zeroes of f2 is exactly 
D + tD. Equivalently, D + tD = 2oo — K, where oo denotes the 
marked point. Regular real solutions are generated by the data with 
the following extra constraint: 

The pair (F, r) is of separating type, i.e. after removing all real ovals 
r splits into 2 components. 

An important example of "solvable" inverse spectral transform is 
the one-energy problem for the two-dimensional Scrodinger operator 
started in the works 13]. 



It is well-known, that the full set of the scattering data for multidimen- 
sional Schrodinger operators > 1 is overdetermined. A lot of people 
studied this problem. We even don't quote this literature. However, the 
case n = 2 turned out to be very specific. In 1976 Manakov, Dubrovin, 
Krichever and Novikov ^Hl El started completely new approach for 
this specific case: They started the Inverse Scattering theory and cor- 
responding Soliton Theory associated with one selected energy level. 
A lot of work was done later in this subject later (see [33 EH UHI and 
review ^J.^ for additional references). In particular, in the first work 
[S] they defined the natural analogs of finite-gap potentials for the 2- 
dimensional problem as the potentials, "finite-gap at one energy" . Let 
u{x, y) be double-periodic. Denote the dispersion relation by ejik^, ky). 
The Fermi-curve at the energy level Eq is defined by: 

(16) €j{k^,ky) = Eo 

Denote the complex continuation of the Fermi-curve by F. The Poten- 
tial u{x, y) is called finite-gap at one energy, if F has finite genus. 

For generic spectral data the operators constructed in |3] have gener- 
ically a non-zero magnetic field, i.e. they have some extra first-order 
terms: 



It might happen that H{x,y) 7^ 0, where H{x,y) = dxA2{x,y) — 
dyAi{x,y). For physical applications it is important to select the case 
of "potential operators" Ai{x,y) = A2{x,y) = with real potential 
u{x,y). Sufficient conditions on the spectral data leading to the po- 
tential operators were found in 1984 by Novikov and Veselov in [S^ . 
For double periodic potentials the existence of such form is necessary. 
It follows from the direct spectral theory, developed by Krichever j^. 



(15) 



L = -dl~ dl + u{x,y) 



(17) 



L = -dl-dl + Ai(x, y)d^ + ^2(0;, y)dy + m(x, y) 
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For the generic regular quasiperiodic potentials "finite-gap for one en- 
ergy level" , this problem remains open. Selection of real potentials here 
is simple. 

How to select the class of regular potentials in terms of algebro- 
geometrical spectral data? There is no complete solution to this Prob- 
lem. It was shown in [33], that if the spectral curve is the so-called 
M-curve, then the potential u{x,y) is regular, and the operator L is 
strictly positive (the selected energy level lies below the ground state). 
An alternative proof of the last statement was obtained by the authors 
in ^B]. The complete characterization of the data generating strictly 
positive operators (with real regular potentials) was "more or less" 
clarified but some special features remain unproved rigorously. 

If the selected energy level is located above the ground state, the 
topology of the spectral curve F become more complicated. Many 
classes of spectral data generating real non-singular solutions were 
found by Natanzon (see review f^), but the classification is not com- 
plete till now. 

3. Sine-Gordon equation. 

Connections between the Sine-Gordon equation and the inverse scat- 
tering method were first established by G.Lamb in 1971 The mod- 
ern approach developed by Ablowitz, Kaup, Newell and Segur in 1974 
[S] is based on the following zero- curvature representation: 

(18) vl>^ = l(f/ + V^)v[/, v[/, = l(f/-V)vl/, 
where 

(19) U = U{X,x,t)-- 



i{ux + ut) 1 
-A -i{ux + Ut) 



(20) V = V{X,x,t) 



-^e™ 
e"^" 0. 



As we mentioned above, the finite-gap "spectral data" consist of 

(1) A hyperelliptic Riemann surface F with 2g + 2 branching points 

(0,Ai,...A2<„oo): /i2 = An(A-A,) 

i=l 

(2) The divisor (a collection of points) D = 71 + . . . + 7^ in F. 

In our text we always assume, that the spectral curve F is generic, 
i.e. all branching points are distinct. 

Construction of the complex Sine- Gordon solutions is based on the 
following standard Lemma: 
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Lemma 1. For generic data T, D there exists a unique two- component 
V ector-f unction ^>{'^,x^t) (the "Baker-Akhiezer" functions) such that 

(1) For fixed {x,t) the function \l/(7,x, t) is meromorphic in the 
variable 7 G F outside the points 0, 00 and has at most 1-st 
order poles at the divisor points •jk, k = 1, . . . , g. 

(2) \l/(7,x,t) has essential singularities at the points 0, 00 with the 
following asymptotic: 

^ ^ V , ; iVXMx,t) +0{\) J 

with some 0i(x, t), (j)2{x,t). 
The Sine-Gordon potential u{x,t) is defined by: 

<p2{,x,t) 



(23) u{x,t)=ihi 



0i(x,t)' 




Denote by \k{x, t) the projections of the zeroes of the first component 
of \l/(7,x, t) to the A-plane. Then 

(24) e^"(-'*) = n(-A,(x,t)) 

i=o 

Remark. To be more precise, the formulas p8H24j) define simul- 
taneously a pair of Sine-Gordon solutions Ui(x,t), U2{x,t), depending 
on the choice of the branch 1/ ^/{X) near the point A = 0. They are 
connected by the following relation U2{x,t) = Ui(t,x) + n. In the real 
case it is possible to fix a canonical branch by making the analytical 
continuation along the real line. This rule is unstable in the following 
sense: if we add a pair of complex conjugate branching points which 
are very close to the positive half-line (or, equivalently, open a resonant 
point), it is a small transformation in terms of the spectral data, but 
it exchanges Ui with U2-- 

The real Sine-Gordon solutions (by Cherednik's lemma they are au- 
tomatically regular |3]) correspond to the following data: 

(1) r is real, i.e. the branching points of F are either real, or form 
complex conjugate pairs. Therefore we have an antiholomorphic 
involutions r : (A, /i) (A, /i). Denote the number of real finite 
branching points by 2 A; + 1. 
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(2) All real branching points lie in the negative half-line A < 0. 
It is convenient to use following enumeration for the branching 
points different from and oo: 0>Ai>A2>...> \2k, 

^2k+l = ^2k+2, • • • 5 -^29-1 = ^2g- 

(3) There exists a meromorphic differential Q (Cherednik differen- 
tial) with first order poles at 0, oo, holomorphic on r\{0, oo} 
with the zeroes at the points 71,. . . ,jg, tji,. . . ,T'yg (or, equiva- 
lently the divisor D satisfy the relation D + tD = + 00 — K). 

As it was shown in |4j , the variety of all real potentials corresponding 
to the given spectral curve F consists of 2^ connected components. A 
characterization of these components in terms of the Abel tori was 
obtained in [7] but this technique did not led to the the calculation of 
topological charge through the inverse spectral data. 

Our calculation of the topological charge for the finite-gap Sine- 
Gordon solutions is based on the following effective description of these 
components (see [12], [EI, [T^l: 

Any meromorphic differential with first-order pole at 00 can be writ- 
ten as: 

^^^^ i^-c^i ^^^^^/^ J 

where Pg_i(A) is a polynomial of degree at most g — l. It is also natural 
to put c = 1. In case of the Cherednik differentials the set of zeroes is 
invariant with respect to r. Therefore all coefficients of the polynomial 
Pg-iW are real. 

Assume, that we take an arbitrary real polynomial Pg_i(A). Is it 
possible to construct a real Sine-Gordon solution corresponding to it? 
The necessary and sufficient condition is the following: the zeroes of 
Q can be divided into two groups {71, . . . , 7^} and {7g+i, . . . , 72^} 
such, that T'-fk = lk+gi k = l,...,g. Equivalently, a polynomial 
Pg-i{X) generates real SG solutions if and only if all real root of Q 
have even multiplicity. In generic situation (all roots form distinct 
complex conjugate pairs) each polynomial Pg-i{X) generates 2^ differ- 
ent solutions. To choose one of them one has to say, which point to 
choose in each complex conjugate pair belonging to D ( the second one 
belongs to tD ). In degenerate cases (i.e. if there are real roots) the 
number of choices is smaller. Al these solutions associated with a given 
Pg-i{\) belong to the same real Abel torus. 

Definition. A polynomial Pg_i(A) (and the corresponding differen- 
tial Q) are called admissible if all real roots of Q have even multiplicity. 

Admissible polynomials Pc,_i(A) can be characterized in the following 
geometrical way: 
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(26) MX) = 

and fill in the following domains by the black color: 



X<0,y'< 



R{X) 



2 5 




s„=+l s =+1 s,=-l 




Fig 1. 

Lemma 2. The polynomial Pg_i(A) is admissible if and only if the 
graph of Pg^i{X) has no parts inside the black open domains. 

If the graph does not touch these domains, we have no real divisor 
points. Real divisor points correspond to the case, when the graph 
touches one of these domains but does not cross the boundary. 

Each pair \2j-i-, is connected by a black "island". The graph of 
admissible Pg_i(A) should go above or below this island, therefore at 
all intervals [X2j, ^2j-i], j < ^ Pg-iW 7^ 0. Let us associate with an 
admissible polynomial Pg_i(A) a collection of numbers sj, j = 1, . . . , k 
by the following rule: sj = 1 if the graph of Pg_i(A) is positive in 
the interval [)^2j, ^2j~i] and Sj = —1 otherwise. Let us call the set Sj 
topological type of the real solution. We have exactly 2^ possible 
topological types. Elementary analytic estimates (see [T^) show, that 
all these components are non-empty. Each connected component is 
a real Abel torus, and the x-dynamics defines a straight line in this 
torus. To calculate the density of the topological charge it is sufficient 
to know the direction of this line and the charges along the basic cycles. 
It follows from a simple analytic lemma: 

Lemma 3. Let u{X), X G M" be a smooth function in R" such, that 
exp(m(X)) is single-valued on the torus ]R"/Z". Equivalently it means, 
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that exp{iu{X + A^)) = exp{iu{X)) for any integer vector N , and 
(28) 

u{X\ X^, . . . , X'^ + 1, . . . , X") - u{X\X\ . . . , X^ . . . , X") = 27mk. 

The numbers are called the topological charges along the basic 
cycles 2tfc, /c = l,...,n. Denote by u{x) restriction of u{X) to the 
strait line X = Xq + x ■ v, v = (t'^,f^, . . . ,f"). Then the density of 



topological charge n 
does not depend on the point Xq and: 



lim [u{x + T) — u{x)\/2'kT is well-defined; it 

T — *oo 



(29) 



n 



fc=i 



The calculation of the direction vector for the x-dynamics is abso- 
lutely standard (see, for example [lOj). Denote by cu' the canonical 
basis of holomorphic differentials on F: 



(30) 



UJ 



9-1 



Then for the components of the x-direction vector we have: 

(31) f/fc = 

To obtain a simple expression for the basic charges it is critical to 
use a proper basis of cycles in F. In [in]-|IS]) the authors used the 
following basis, first suggested in P : 





bk=[^2k.^2k-llk=l,2,3 



Fig 2. 
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Here the cycles aj, j = 1, . . . , are ovals on the upper sheet of F, 
containing inside the points Aq = 0, Ai, A2,. . . ^2j-i- The cycle bj, 1 < 
j < k lies over the interval [\2j, A2j-i]- The cycles Qj, j = k+1, . . . ,g\ie 
over pathes connecting the pairs X2j-i and X2j- We assume that these 
cycles do not intersect each other, and the cycles aj, j = k + 1, . . . , g do 
not intersect the negative semi-line. The cuts are shown by the zigzag 
lines. The upper sheet contains the semi-line A > 0,/i > 0. 

Consider a basic cycle 2tj on the real component of Jacoby torus, 
represented by the closed curve. The image of this cycle in F under 
the inverse Abel map is a closed oriented curve Cj, formed by the mo- 
tion of the corresponding divisor points (it may have several connected 
components). The motion of an individual divisor point does not have 
to be periodic, after going along the cycle we may obtain a permu- 
tation of the divisor points. The curve Cj is homological to the cycle 
Qj G i^i(F, Z). It follows from that the topological charge rij along 
the cycle 21^ equals to the winding number of the curve Cj with respect 
to the point 0. Equivalently 



(32) nj = CjoR_, 



where o denote the intersection number, Cj denotes the projection of Cj 
to the A-plane, R_ is negative semi-line with the standard orientation. 

For each point of % the corresponding divisor 71,. . . ,7^, is admissible. 
From the characterization of admissible divisors obtained above it is 
easy to show, that the curve Cj does not touch the closed segments 
on the real line [—oc, \2m],- ■ ■ ,[^3, A2], [Ai, 0]. Therefore any time the 
curve Cj crosses the negative semi-line, it intersects one of the basic 
cycles bj, j = l,...,k. 

Unfortunately this information is not sufficient to calculate the ba- 
sic charge, because the orientation of the cycles bj coincides with the 
orientation of the negative semi-line at one sheet and they are oppo- 
site at the other one. For example at the Fig. 3 we see two different 
realizations of the cycle ai representng different topological types, ai 
is drawn at the upper sheet and a'l is drawn at the lower one. We have 
ai o bi = a'l obi = 1, but ai o R_ = 1, d[ o M_ = —1, therefore ni = 1 
and ni = —1 for these cycles respectively. 
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a! 



^2 




Fig 3. 



Fortunately, the topological type contains the information at which 
sheet the intersection takes place. Namely we have: 

Lemma 4. Assume, that the cycles Cj intersects the negative semi-line 
at the interval {X21, ^21-1) ■ Then orientations of bi and M_ coincide in 
the intersection point if {—lY^^Si > and are opposite if {—iy~^Si < 0. 

Combining all these results we obtain the final formula: 

Theorem 1. The density of the topological charge for a real Sine- 
Gordon solution is given by 

k 
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